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In this paper we investigate the influence of low frequency turbulence on Doppler spectral line 
shapes in magnetized plasmas. Low frequency refers here to fluctuations whose typical time scale 
is much larger than those characterizing the atomic processes, such as radiative decay, collisions 
and charge exchange. This ordering is in particular relevant for drift wave turbulence, ubiquitous 
in edge plasmas of fusion devices. Turbulent fluctuations are found to affect line shapes through 
both the spatial and time averages introduced by the measurement process. The profile is expressed 
in terms of the fluid fields describing the plasma. Assuming the spectrometer acquisition time to 
be much larger than the turbulent time scale, an ordering generally fulfilled in experiments, allows 
to develop a statistical formalism. We proceed by investigating the effects of density, fluid velocity 
and temperature fluctuations alone on the Doppler proflle of a spectral line emitted by a charge 
exchange population of neutrals. Line wings are found to be affected by ion temperature or fluid 
velocity fluctuations, and can in some cases exhibit a power-law behavior. This study gives some 
insights in the appearance of non-Boltzmann statistics, such as Levy statistics, when dealing with 
averaged experimental data. 
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I. INTRODUCTION 

Spectral line shape studies have played a major role 
in the investigation of the nature of atomic radiators 
and their environment, in astrophysics as well as in lab- 
oratory plasmas. Indeed, depending on the dominant 
line broadening mechanisms, it is for instance possible 
to retrieve the electron density or the ion temperature 
from the analysis of a given line. However, in many 
cases actual plasmas are far from thermal equilibrium, 
being inhomogeneous or having non-Maxwellian velocity 
distributions, features which significantly complicate the 
analysis of experimental data. In addition, these depar- 
tures from thermal equilibrium can trigger instabilities, 
whose growth and non-linear saturation eventually lead 
to the onset of turbulence J* i2„ ^] . The importance of 
investigating the possible effects that turbulence might 
have on line shapes has been acknowledged very early, 
and the motivations of these studies were, and still are, 
two-fold : first there is the need to quantify the errors 
introduced by neglecting turbulence in routine diagnosis 
based on line-shapes. Then, the possible existence of sig- 
nificant deviations could be used to diagnose turbulence 
itself. An optical diagnostic of turbulence based on pas- 
sive spectroscopy would indeed be very convenient, this 
sustainable technique being already available on numer- 
ous experiments. 

Historically, a large number of papers (0,0 and ref- 
erences therein) have dealt with the Stark effect result- 
ing from turbulent electric fields, such as those associ- 
ated to a high supra-thermal level of Langmuir waves Q . 
Starting from the seminal paper by Mozer and Baranger 
0, several models have been devised to include turbu- 
lent Stark broadening in the calculation of line-shapes. 
The results thus obtained are relevant to plasmas for 



which Stark effect is dominant compared to Doppler ef- 
fect. There are however situations for which this ordering 
is reversed, and important examples are edge plasmas of 
magnetic fusion devices such as Tokamaks in the ionizing 
regime. For these low density plasmas (A^e < 5 x 10^° 
m~'^), Zeeman and Doppler effects are the dominant 
broadening mechanisms for low lying lines such as the 
Da (transition between the levels n = 3 and n = 2 of 
the atomic deuterium). In such cases, line shape stud- 
ies essentially provide measurement of the emitters ve- 
locity distribution, and have so far brought valuable re- 
sults concerning the origin of neutrals in edge plasmas 
d, 0, 0, 01 . However, these plasmas are known to be 
strongly turbulent, i.e. the level of fluctuation of the fluid 
fields characterizingthe plasma can rise up to several 
tenths of percents [J- The experimental characteriza- 
tion of these fluctuations is of first importance to analyse 
drift-wave (DW) turbulence, which is held responsible for 
the so called anomalous transport degrading the quality 
of the confinement Q . 

As an example, we will consider the case of the Balmer 
a of hydrogen isotopes (Da for the case of deuterium), 
since it is one of the most routinely monitored line in 
edge plasmas, being both intense and optically thin. In 
sections ^ and IIIII the expression of the measured line 
profile is carefully discussed to emphasize the role of the 
spatial and time averages involved in the measure. We 
show in section ^] that a neutral population created by 
charge exchange can be considered as being in a local 
equilibrium characterized by the local density, tempera- 
ture and fluid velocity of the ions. In sections IIVI and 
fVl we will show that in presence of low frequency turbu- 
lence, the Doppler profile gives access to an apparent ve- 
locity distribution. By further developing the model only 
briefly presented in [3j j ^^^^ apparent VDF is reex- 
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pressed in terms of the Probability Distribution Function 
(PDF) of the fluid fields. In section IVIl the influence of 
density, fluid velocity and temperature are successively 
investigated in details. Finally, it is shown in section IVlTI 
that for particular choices of the statistical properties of 
the turbulent fluctuations, the apparent VDF becomes 
a Levy distribution. This result establishes a clear con- 
nection with one of our previous work [14] . in which we 
investigated the possible origin of a power law behavior 
observed in the line wings of Da spectra measured in the 
former crgodic divertor configuration of the Tore Supra 
Tokamak. 



II. EXPRESSION OF THE MEASURED 
SPECTRA 

Let us first define precisely the observable quantity 
for a spectrally resolved passive spectroscopy measure- 
ment. First of all, obtaining the spectrum emitted by 
the plasma (which will be referred to as the measured 
spectrum in the following) from the raw spectrum in- 
volves deconvolution of the apparatus function Tap. In 
practice, the theoretical spectrum is convolved with Tap 
before being compared to the raw spectrum. The ra- 
diation emitted by the plasma is integrated both along 
the Line Of Sight (LOS) and during the acquisition time 
of the spectrometer, denoted by r^. The observable in- 
tensity X,„es(AA), where AA stands for the wavelength 
detuning from the center of the line, is thus given by the 
following expression 

J„e.(AA)-/ TiociAX,z,t)—^dz dt, (1) 

Jo JC 4770 

where X/oc(AA, z,i) is the local line shape emitted at a 
given distance z from the detector along the LOS £. Here 
6S stands for the detector active area. Assuming the 
latter to be delimited by zi < zi such that L — zi — z\ <SL 
zi, Eq. J^l reduces to 

X™es(AA)~— ^/ / TiociAX,z,t)dz dt. (2) 
We now introduce the local absolute brightness 

/-|-cx> 
Tioc{AX,z,t) dAX, (3) 
-oo 

and define the local line shape normalized to unity by 
T{AX, z,t) — Tioc{AX, z,t)/b{z,t). In the remainder of 
this paper, we will deal with the measured profile nor- 
malized to unity, given by 



Tmes{AX) — 




B{z,t)T{AX,z,t)dz dt. 



(4) 



where the relative brightness B is defined by 

— Jo Ljz,^i^^*)dtdz 

So, the measurement process both entails a spatial 
and a time average of the local profile. In order to 
achieve time resolved measurements, the acquisition time 
Tm should be chosen shorter than the typical turbulent 
time, which is associated to the time variations of the 
functions B{z, t) and T{AX, z, t). However, such a choice 
would generally result in spectra having very low signal to 
noise ratios. To retrieve information from time resolved 
measurements, one must then either forsake spectral res- 
olution or resort to active techniques, such as Beam Emis- 
sion Spectroscopy (BES), which allows to diagnose the 
time behavior of turbulent density fluctuations at the 
edge of tokamaks (e.g. HE El)- Conversely, if spectral 
resolution is needed, the acquisition time should be suf- 
ficiently large so as to ensure reasonable signal to noise 
ratios. In this paper we will deal with the latter situa- 
tion, for the limiting case in which the acquisition time 
is much larger than the typical turbulent time scale. In 
fact, this generally corresponds to the actual situation for 
passive spectral line shapes measurements. Spectra thus 
obtained will prove to yield further information than time 
resolved experiments. The next step in the modelling 
consists in relating the local brightness B and the local 
profile T to the parameters characterizing the plasma. 

III. MODELLING OF THE LOCAL PROFILE 

In this section we will first describe the model that 
will be used to describe the plasma, i.e. a set of fluid 
equations. The remainder of the section will present the 
expressions of the local brightness and profile relevant to 
edge plasmas typical conditions. 

A. Plasma description 

We are interested in plasmas which can be described 
by a set of N macroscopic fields denoted by X = 
{Xi(r, t), . . . , XAr(r, t)}, including the density, temper- 
ature and fiuid velocity for each species. These fields are 
solutions of a set of fluid equations. For example, the 
field Xi{v,t) would satisfy to a conservation equation 

BX 

^ + v.r = 5(x), (6) 

where ^(X) is a source term and F is a flux, for exam- 
ple given by the sum of a diffusive and a convective flux 
r = —nVXi + uXi, u = Xj being a velocity field. It 
is furthermore assumed that this set of fluid equations 
describes a turbulent stationary state, for which the sta- 
tistical properties of the plasma do not change during the 
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acquisition time Tm- For each species, the vahdity of a 
fluid description rehes on the ordering r ^ z/^J^ , where r 
is the typical time of variation of the fluid fields and Vcoii 
the collision frequency (in the case of DW turbulence in 
edge plasmas, we have t ~ 10 — 100 /is). The fulfillment 
of this ordering ensures that the VDF of the species 
7 = i, e remains close to a local Maxwellian at each time 
and location 



2nT^{r,t) 



cxp 



mj{-v — u-y(r, t))'' 



(7) 

where m^, Ty(r, t) and u^(r, t) are respectively the mass, 
the temperature field expressed in eV, and the fluid ve- 
locity fleld of the species labeled by 7. In the following, 
we shall consider the case of a pure deuterium plasma for 
which Zeff = 1. 

The calculation of the neutrals VDF requires the use 
of a reflned model. Indeed, there are different sources of 
neutrals in edge plasmas of tokamaks, each of them giving 
birth to a single class of neutrals. These classes, charac- 
terized by different temperatures, coexist since the den- 
sity is usually too low in order to ensure their complete 
relaxation toward the background local equilibrium. The 
lowest temperature class originates from the dissociation 
of molecules released from the wall, whereas those having 
larger temperatures are mainly attributed to charge ex- 
change reactions (e.g. [1, IE Cfl Il3 ) ■ ^he following, we 
will only consider the class of neutrals locally created by 
charge exchange reactions, that plays an important role 
for the line wings behavior. Indeed, it will be shown that 
turbulence essentially affects these regions of the spec- 
tra. In order to model the VDF of these emitters, we 
can once again take advantage of the separation of scales 
between atomic processes and turbulence. In fact, the 
inverse of the charge exchange rate is of the order of a 
few /xs, i.e. shorter than the typical turbulent time scale. 
As a result, the emitters VDF remains at each time close 
to that of the ions, given by Eq. O with 7 = 1. From 
the microscopic point of view, the emitter's VDF thus 
appears as a Maxwellian characterized by a set of slowly 
varying macroscopic flelds. 



B. The local brightness 

The local line brightness is directly related to the 
population of the transition upper atomic level. In 
general, this population has to be calculated by tak- 
ing into account the contributions of the different pro- 
cesses (for instance collisions, charge exchange, radia- 
tive decay) populating or depopulating the levels. If 
the fluid fields characterizing the plasma vary slowly on 
the typical time scales associated to these processes, a 
stationary approach is suitable to calculate the bright- 
ness. The levels populations are assumed to be time 
independent and are calculated using the values of the 
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FIG. 1: plot of the brightness per emitter as a function of 
the electronic temperature Te in the range to 100 eV for a 
given density Ne = 10^^ m~^. The temperature dependence 
is weak for Te > 15 eV. 



fields X(r,i) at each time and location. In practice, 
the brightness essentially depends on the electron density 
Ne{r,t) and temperature Te(r,t). We have performed a 
calculation of the brightness per emitter Bi (defined as 
B{r, t) = no(r, t)Bi(K{r, t)), where no(r, t) is the density 
of emitters) for the Da line in edge plasma conditions, 
i.e. Ne = lO^s - 10^3 m~3 ^^^j = 1 - 100 eV, using 
the code SOPHIA jl3|. The electron density dependence 
of Bi is found to be linear, in accordance with the fact 
that the upper level of the transition is essentially popu- 
lated by electronic collisions from the ground state. This 
leads to a quadratic behavior of the brightness with Ne, 
since no Ne- The infiuence of the electron tempera- 
ture on the brightness is more subtle, as shown on Fig. 
n for two different densities. The existence of a maxi- 
mum reflects the competition between the growth of the 
electron collisions cross section with temperature, which 
dominates the small temperatures behavior, and the ion- 
isation process. For electron temperatures larger than 15 
eV, the influence of Te on the brightness is weak, and in 
the remainder of this paper, we shall therefore consider 
the brightness as being only a function of the electron 
density. 



C. The local line shape 

The AA dependance of the local profile T{AX,z,t) 
is determined by the dominant line broadening mech- 
anisms. In magnetized plasmas, Zeeman, Stark and 
Doppler broadenings should a priori be taken simultane- 
ously into account. In general, the local profile normal- 
ized to unity can be written as the following convolution 
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product 



I{AX,z,t)= J dAX'IzsiAX- AX',z,t)lD{A\,z,t), 

(8) 

where Izs is the local Zeeman-Stark profile, which de- 
scribe the broadening resulting from the effect of the 
magnetic and electric fields on the emitters energy levels 
P|. The Doppler profile Id is related to the wavelength 
shift introduced by the movement of the radiator along 
the LOS, and is thus directly given by 

lD{AX,z,t)dAX^ f{v„z,t)dv„ (9) 

where f{vz,z,t) stands for the emitters VDF along the 
LOS, obtained from |2J upon integrating over the two 
components of the velocity perpendicular to the LOS 



f{vz,z,t) = J J dvxdvyF{-v,z,t). (10) 

It should be noted that Eq. ^ would not be valid if 
the velocity of the emitter were not constant during the 
emission process, due to collisions Ts'l. If Aljo denotes 
the Doppler line width expressed in units of pulsation, 
Eq. © assumes that t^J^^ ^ Auj^. This ordering is 
largely satisfied in edge plasmas, and is moreover not in- 
consistent with the assumption TcoI > t underlying the 
validity of Eq. Q. For a given line, the relative impor- 
tance of the different broadening mechanisms depends 
on plasma conditions, i.e. on the average values taken 
by the plasma density and temperature, but also on the 
detuning AA. In the following, we will again discuss the 
case of the Da line, first for the bulk of the line and 



where m is the emitters mass, Aq the unperturbed wave- 
length of the transition under study, T{r, t) the ion tem- 
perature, and Uz(r,t) the component of the ion fluid ve- 
locity along the LOS. 



IV. APPARENT VELOCITY DISTRIBUTION 

Gathering the results of the above sections, we obtain 
the following expression for the measured profile normal- 



then for line wings, these regions of the spectra for which 
|AA| > AAi/2, AAi/2 being the HWHM of the profile. 
In the center of the line, Stark effect is negligible for den- 
sities lower than TVg = 5 x 10^" m~'^, an ordering which 
is usually (but not always) satisfied in edge plasmas. In 
addition, since the magnetic field is larger than 1 T, fine 
structure can be neglected fig^. Therefore, the Dq line 
splits into three Doppler-broadened Zeeman components 
(one TT and two a) . The lateral a components are equally 
separated from the central tt component. Under parallel 
observation with respect to the magnetic field, only the 
cr components are observable. Although negligible in the 
bulk of the line. Stark effect might become dominant in 
the line wings for detunings larger than a value AXs{Ne) 
which is an increasing function of the density. Therefore, 
in the remainder of the paper it should be understood 
that the Doppler line wings are the regions of the spectra 
for which both orderings |AA| ^ AA1/2 and |AA| < AA5 
are simultaneously valid. The existence of such a regime 
depends on the plasma conditions. Its study is relevant 
for edge plasmas and consequently Stark effect will be ne- 
glected in the remainder of the paper. However, it should 
be emphasized that the statistical formalism which is de- 
veloped in section would also be applicable if Stark 
effect were not negligible. In the latter case, the local 
profile should be calculated using Eq. JSJ instead of Eq. 

According to Eq. Q , the Doppler spectrum of a single 
Zeeman component is proportional to the emitters VDF 
/ along the line of sight. As previously explained, we con- 
sider a class of neutrals created by charge exchange reac- 
tions, whose VDF is approximated by a local Maxwellian. 
The corresponding expression of the local Doppler profile 
is given by 



(11) 

I 

ized to unity 

I^es(AX) ^ — f" dty f dzB{X{z,t))lD{AX,X{z,t)), 

Tm Jo ^ J C 

(12) 

which is now expressed in terms of the fluids fields de- 
scribing the plasma. The apparent velocity distribution 
function fa(vz) is straightforwardly deduced from the 
measured spectrum by 

I,nes{AX)dAX = .fa{v,)dv,, (13) 

in analogy with Eq.@. This VDF is an average of the 



I 

In{AX,Tir,t),u4r,t)) = ^^^^exp ^^^^ j , 
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local emitters VDF over time and space. Indeed, combin- 
ing Eq. (|12|l and Eq. leads to the following explicit 
expression 

/aK) = — r^dty [ dz B{X{z,t)) f{v,,X{z,t)). 
Ttti Jo ^ J c 

(14) 

The apparent VDF /„ can be given a deep physical 
meaning as will be shown in section FVIII 

Intuitively, in plasmas where the fluctuation rate is low, 
fa should remain close to a Maxwcllian f^q characterized 
by the time and space averaged values of the temperature 
and the velocity fields, respectively denoted by T and u^, 
i.e 

fa{Vz) - feq{Vz;T,U^). (15) 

Conversely, in a situation where strong fluctuations oc- 
cur, there is a priori no obvious reason for which the 
apparent velocity distribution should remain close to the 
average Maxwellian given by Eq. (|15|l . In particular, 
in edge plasmas the fluctuation rate can rise up to sev- 
eral tenths of percents. The validity of Eq. IjlSI) then 
clearly becomes questionable, and Eq. H14I) should be 
used instead. A calculation of the apparent VDF fa can 
be carried out from the latter equation once the solu- 
tions of the fluid equations are known, i.e. the time and 
space dependences of each of the fields Xi{r, t) have been 
worked out. Due to the non- linear nature of the fluid 



equations and the complexity of the geometry, this cal- 
culation would best be achieved numerically. Although 
such an approach might be able to encompass the com- 
plexity of the problem, we find it worthwhile to begin 
with a simpler one in order to gain insights on the kind 
of effects that turbulence might produce on spectral line 
shapes. 



V. STATISTICAL FORMALISM 

A. Expression of the profile 

In the following, we will take advantage of the fact that 
the acquisition time of the spectrometer is usually much 
larger than the typical time scale of the turbulence r. 
Let us first note that upon using an appropriate normal- 
isation for the 5 function, the following relation holds for 
any z and t 



[ f[S{X,-X,{z,t))dXi..dXN^l, (16) 

where Xi is the sample space variable corresponding to 
the Held Xi{z,t). Introducing this identity into Eq. H14|l . 
interchanging the order of time and sample space integra- 
tions, and finally making use of the delta function sifting 
property yields the following expression for the apparent 
velocity distribution function 



fa{Vz) = T I '^^ I dXi...dXN 

N 



N 



— \{5{X,-X,{z,t))dt 



B{X) fivz,-K). 



(17) 



The quantity between brackets is a time average of the delta functions product, whose typical time variations occur 
on the time scale t ^ r^. It is therefore justified to let tend to infinity |23|, and then use the ergodic assumption, 
i.e. replace the time average by an ensemble average denoted by the brackets (•) 



lim — / \\5{X,-X,{z,t))dt=lT\5{X,-X,{z,t)) 



(18) 



i=l 



\i=l 



r 



This ensemble average has to be understood as an av- 
erage over the time realisations of the stochastic fields 
Xi{t,z) at point 2, assumed to be a stationary process. 
Introducing the local joint Probability Density Function 
(PDF) of the fluctuating fields defined by 



N 



V{Xi,--- ,Xm,z) = [\{5{X,-X,{z,t)) 



(19) 



the apparent VDF becomes 



fa(Vz) = - dz dXViX,z)B{X) fivz,X). (20) 
^ Jc Jn 



Finally, upon integrating on the space coordinate z, the 
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apparent VDF is given by 



fa{v,)= f dXW{X)B{X) f{v,,X), (21) 

where the spatiahy integrated PDF VF(X) is obtained 
from 



W{X) ^-J dzV(X,z). 



(22) 



In the remainder of the paper we shaU furthermore as- 
sume homogeneous turbulence, that is M^(X) = P(X, z) 
(note that the weaker assumption of homogeneity along 
the line of sight is sufficient). 



B. Discussion 

In the frame of our statistical reformulation, it is no 
longer necessary to know the solutions of the fluid equa- 
tions in order to calculate the apparent VDF. Instead, 
the joint PDF of the turbulent fields should have been 
computed. A straightforward approach would be to rely 
on a fluid code, so as to compute histories of the dif- 
ferent fields, and then their PDF. As we have already 
pointed out, this would require heavy numerical compu- 
tation, especially in order to obtain the PDF tails with 
a good accuracy. Furthermore, if such calculations were 
carried out, any statistical reformulation would obviously 
be superfluous, and the apparent VDF could directly be 
obtained from Eq. An approach more suited to 

our formalism should proceed directly at the PDF level. 
The next section will be devoted to present such a model, 
initially developed by Pope However, it should be 

emphasized that such a calculation is bypassed if assump- 
tions for the shape of the PDF are made. This is the one 
of the advantages of our formalism, since it allows to draw 
conclusions on the properties that turbulence should have 
so as to significantly affect line shapes. 



C. Determination of the PDF from the fluid 
equations 

Let us consider the passive advection of a scalar X{z,t) 
solution of Eq. lO, in which the source term is an ar- 
bitrary function of X and the flux T is the sum of a 
convective term and a diffusive term. The convective ve- 
locity field u is assumed to be an incompressible stochas- 
tic field, the statistical properties of which are known. In 
order to calculate the apparent velocity distribution func- 
tion from Eq. H21|l . the spatially integrated joint PDF of 
velocity and temperature, denoted by W{u,X), should 
be calculated. Here, we will limit ourselves to the mod- 
elling of the marginal distribution W{X), obtained by 
integrating W{u,X) over the velocity. Indeed, this will 



be sufficient to highlight the salient points of the model. 
Assu ming homogeneous turbulence, and then following 
Pope |2l| . the time dependent PDF W{X,t) is shown to 
obey a Fokker-Planck like equation 



dW d - 5^ - 

(23, 

where S{X) is the source term in the fluid equation. The 
expression of the function D{X) will be discussed below. 
The stationary solution of the latter equation is 



W{X) 



c 



D{X) 



exp 



X 



Sjw) 
D{w) 



dw 



(24) 



As a result, in the PDF approach a non- linear source 
term S does not introduce any closure problem, unlike 
in the moment based models [l|. The problem remains 
nonetheless unclosed, since the function D{X) is in gen- 
eral not expressible in terms of W{X) or S{X) alone. 
Indeed, the shape of this function depends on the cor- 
relations between X and its gradient. More precisely, it 
can be recast in the following form 



D{X) = 



d{WX)V{WX\X){WXf, (25) 



where ViVXlX) is the PDF of the gradient of X, condi- 
tioned to a given value of X 21]. In order to obtain this 
PDF, an equation for the joint PDF of X and its gradient 
should be written |22l |. which in turn would involve cor- 
relations with higher orders gradients. Eventually, one 
ends up with an inflnite hierarchy of equations, involv- 
ing the joint PDFs of X,VX,V^X, ... . In addition, it 
should be kept in mind that the statistical properties 
of the velocity field u, while not appearing explicitly in 
Eq. 123), do actually affect the shape of D{X) through 
Eq. as should the expression of S{X). The closure 
of this hierarchy has proven to be difficult to address. 
Promising techniques, such as the mapping closure 
have been devised to overcome these difficulties, but have 
not yet led to decisive results (for an application to the 
Hasegawa-Mima equation governing plasma turbulence, 
see Ref. jl^). Addressing these issues is largely beyond 
the scope of the present paper, and for our purposes 
it will be sufficient to present an early attempt to this 
closure problem, due to Sinai and Yakhot h4. These 
authors were interested in the case of passive advection 
of temperature in homogeneous decaying turbulence, for 
which there is no source term in the temperature equa- 
tion. Their idea is to deal with the rescaled quantity 
X — T I (T^), which is solution of an equation analogous 
to (jHJ, S{X) being a linear function of X. The following 
Taylor development is used to express the function D{X) 



D{X) ~ l + fcV^ 



(26) 
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FIG. 2: Plot of the Sinai PDF for To = 30 eV, it = 10 and 
for different values of ct = 5, 10, 20 eV. These distributions are 
used to compute the corresponding apparent VDF on Fig. |^ 



VI. APPLICATION TO THE CASE OF ONE 
FLUCTUATING VARIABLE 

In an actual turbulent plasma, several fields fluctuate, 
and these fluctuations are coupled. According to Eq. 
(|21|l . the joint PDF of the relevant fields should be com- 
puted before calculating the apparent VDF. However, the 
role of density, velocity and temperature fluctuations on 
the apparent VDF shape have no reason to be identical. 
As a first approximation, it is therefore rational to con- 
sider the idealized case in which only one field fluctuates. 
This will shed light on which field fluctuations lead to 
the most significant effects on line shapes. 



A. Density fluctuations 

Let us first consider density fluctuations. Since the lo- 
cal VDF normalized to unity does not depend on density, 
the integration over density fluctuations is trivially per- 
formed, and the apparent VDF is found to be equal to 
the local emitters VDF 



where the parameter fc > is a measure of the correla- 
tions strength. From Eq. H24|) . the following result for 
the temperature PDF is readily obtained 



W{T) 



C 



l + l/2fc ' 



(27) 



where C is a normalization constant, and a controls the 
width of the distribution. Fig. El shows a plot of W{T) 
for To = 30 eV, = 10 and ct = 5, 10, 20 eV. A few sub- 
tleties and limitations concerning the use of this result 
deserve to be mentioned. First, it should be noted that 
a is actually time dependant. We shall assume here that 
the acquisition time Tm is chosen such that a/ a <C Tm- 
This requires a separation of time scales between the tur- 
bulent fluctuations and the decay of the average quanti- 
ties. Secondly, the correlations are treated using the de- 
velopment given by Eq. (|26l) , which is not valid for large 
values of the temperatures. Our results concerning line 
wings should thus be limited to not too large detuning 
AA. Finally, it should be pointed out that the distribu- 
tion given by Eq. H27|l is a Tsallis distribution (e.g. psj ') 
with q = {1 + 2k)^^. Therefore, in this model, tempera- 
ture fluctuations obey to Tsallis non-extensive statistical 
mechanics 2^ when correlations exist, and to Boltzmann 
statistics for vanishing correlations. 

In the next section, we shall use these results as an 
input for apparent VDF calculations. 



faiVz) = 



+ 00 



dnB{n)W{n)f{v,,T) = f{v,,T). (28) 



Therefore, at this level of approximation, Doppler line 
shapes are not sensitive to density fluctuations. The ap- 
parent VDF should thus remain Gaussian with the tem- 
perature T, whatever the shape of W{n). This is in 
sharp contrast with line brightness time resolved mea- 
surements, which essentially provide information on den- 
sity fluctuations. However, it should be noted that for 
cases in which Stark effect is not negligible, Eq. (|28|) no 
longer holds, since the local line shape then strongly de- 
pends on the density. As we have already pointed out, 
the formalism presented here could nevertheless be used 
upon replacing the local Doppler profile by the total pro- 
file given by Eq. © . 



B. Fluid velocity fluctuations 

Let us now investigate the case in which only the fluid 
velocity fluctuates. In the following, W{uz) stands for 
the PDF of the fluid velocity component along the line 
of sight, and cr^ for its variance. Starting from Eq. H21|l. 
the apparent VDF reduces to 



(29) 



which is the convolution product of W and the local 
Maxwellian. Eq. ^ is a well known result in plasma 
spectroscopy, which is mentioned in classical textbooks 
[23. A shape- independent definition of the apparent 
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temperature T^ff from the profile should proceed from 
its second moment 



)v1dv. 



(30) 



In the fluctuations-free case, the actual temperature of 
the emitters T is recovered, whereas if fluctuations do 
occur the apparent temperature is given by 



T„=T 



UZ 



(31) 



where vth is the thermal velocity corresponding to the 
temperature T. The apparent temperature obtained 
from the Doppler line width is thus not rigorously equal 
to the actual temperature of the emitters. This result 
has already been mentioned by several authors, and was 
actually used in the first models retaining the effect of 
turbulence on Doppler fine shapes [S^- In order to ob- 
tain a ten percents discrepancy between Ta and T for 
deuterium emitters, the fluctuation rate should be of the 
order of thirty percents (i.e. fT„ ~ 0.3 Vth)- This effect 
would be stronger for heavy emitters, since their thermal 
velocity is smaller . This estimation suggests that the 
Da line width is not strongly modifled by fluid velocity 
fluctuations. However, considering only the line width is 
not sufficient. In fact, the line shape, i.e. the apparent 
VDF, is often found to be non-Gaussian, and therefore 
from Eq. so should be the PDF W{u^). Recent 

findings in astrophysical spectra |33| , as well as in 

tokamak plasmas for radial velocity fluctuations [sl] in- 
dicate strong deviations from the Maxwellian, especially 
for line wings. As an illustration, let us consider PDFs 
which have power law tails characterized by an exponent 
a such that < a < 2. It is easily shown from Eq. 
(|29|l that the resulting apparent VDF features a similar 
asymptotical dependence 



faiVz) (X 



1 



|Q-|-1 ■ 



(32) 



Such a behavior would for instance arise if the velocity 
PDF were a Levy distribution (sec section^). Examples 
of these distributions are plotted on Fig. O for a — 
0.5, 1, 1.5 and c" — Vth/W. The resulting apparent 
velocity are plotted on Fig. 01 and exhibit a power-law 
behavior in their tails. However, it should be pointed 
out that for this effect to be observable, large amplitude 
velocity fluctuations of the order of a few thermal velocity 
Vth should actually occur. Finally, it should be noted that 
in a magnetized plasma, the physics underlying parallel 
and perpendicular velocity fluctuations are different. The 
latter are related to the electric field fluctuations through 



3 




0.25 



FIG. 3: Plot of the Levy velocity PDF for a = 0.5, 1, 1.5. 
The fluid velocity is plotted in units of the thermal velocity 
Vth- The parameter c characterizing the width of the distri- 
bution (see section fVl is defined by c™ = Vth/W. 
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FIG. 4: Plot of the apparent velocity distribution correspond- 
ing to Levy velocity PDF in a logarithmic scale. The dotted 
line corresponds to the Maxwellian which would be observed 
in the absence of fluctuations. The velocity is plotted in units 
of the thermal velocity vth- The existence of an algebraic de- 
cay of exponent —a + 1 in the tail of the apparent VDF is 
clearly seen. 



whereas the former can arise from a Kelvin-Helmholtz 
like instability associated to the existence of a perpen- 
dicular gradient of parallel velocity. Changing the orien- 
tation of the line of sight would allow to investigate each 
of these different cases. 



C. Temperature fluctuations 



E X B 

—32—' 



Finally, we consider the case where only the ion tem- 

(33) 

perature fluctuates, and for which the apparent VDF 
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reads 

r + co 

faiv.)^ / WiT) fiv,,T)dT. (34) 
Jo 

The latter is not a convolution product, in opposition 
to the case of velocity fluctuations. To begin with, the 
apparent temperature defined by Eq. H30|l . is given by 

r+oo 

Ta= dT WiT) T, (35) 

Jo 

and is thus equal to the mean temperature of the distri- 
bution W{T). Hence, Tq does not depend on the temper- 
ature fluctuations variance. The profile is obtained as a 
weighted sum of Gaussians of different widths, and thus 
cannot stay rigorously Gaussian itself. Nevertheless, for 
a sharp temperature PDF peaked around Tq, the actual 
deviations from Gaussianity should not be very impor- 
tant, as the dominant contribution in the integral of Eq. 
(|34|l is expected to come from the neighborhood of Tq. 
However, while leading to accurate results for the central 
part of the profile, this line of argument is not correct for 
the line wings. Indeed, the value of /(To, Vz) scales with 
Vz as 

f{To,Vz)(xe^p(^-^y (36) 

and therefore strongly decreases as Vz increases. Con- 
sequently, as shown on Fig. |S1 the contribution of the 
maximum of the temperature PDF in the integral be- 
comes negligible for large enough Vz (i.e. in the wings 
of the apparent VDF), and this especially if W{T) has a 
slowly decreasing tail. 

For instance, an algebraic behavior for the tempera- 
ture PDF implies a similar one for the measured profile. 
The relation between the exponents can be obtained in 
the following manner, noting that for large velocities the 
apparent VDF can be approximated by 

faiv.)- r°^W(T)^dT. (37) 

Using then a power-law ansatz for the temperature PDF, 
the following result is readily obtained 

W{T) « ^ ^ fa{Vz) OC j-^. (38) 

For example, let us consider the case in which the tem- 
perature fluctuations PDF is the Sinai distribution given 
by Eq. 127|l . and plotted in a logarithmic scale on Fig. 
12 for k ~ 10, To — 30 eV, and for different values of a 
ranging from 5 to 20 eV. The bulk of the apparent VDF 
remains very close to that of the Maxwellian at 30 eV 



0.3 




20 40 60 80 100 
T(eV) 

FIG. 5: Plot of the local VDF /(T, v^,) as a function of T for 
three different values of the component of the velocity along 
the LOS Vz = Vth, "ivth, 3vth, where Vth is the thermal ve- 
locity for 30 eV. A model distribution W{T), peaked around 
To — 30 eV is also plotted (solid line). As is increased, the 
contribution of To in the calculation of the apparent VDF be- 
comes all the more negligible than the tail of W{T) decreases 
slowly. 

(dotted line) for every value of a. However, the discrep- 
ancies become important in the apparent VDF tails (i.e. 
Doppler line wings), all the more so a is increased. In 
addition, the tails are found to exhibit a linear behavior 
in logarithmic scale, which signals a power-law depen- 
dence. The exponent which characterizes this algebraic 
decay should take the value —3 — 2/k according to Eq. 
(j2Hl- The k dependence can be checked on Fig. [3 where 
the apparent VDF is plotted for cr = 10 eV and for differ- 
ent values oi k (k = 1, 1/2, 10), i.e. different correlation 
strengths. The stronger the correlations, the larger the 
deviations from the Maxwellian. 

In the frame of the Sinai model, the exponent a char- 
acterizing the apparent VDF power law decay is larger 
than 3. Other turbulence models could lead to smaller 
exponents. Let us indeed investigate the case in which 
the temperature PDF is a Levy distribution of indexes 
< a < 2 and -1 < /3 < 1, denoted by C^A^) HI. In 
the Fourier space, one has 

lnZ„,0(fc) = -clA^r [l + a)) : (39) 

where c controls the width of the distribution, and the 
function Lo{k, a) is defined by 

, , jtan(7ra/2) for a ^ 1, 
"('='")^\(2A)ln|fc| fora = L 

For < a < 1 and /3 = -1, W{T) = for negative 
arguments, as should be the case for the temperature 
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V/Vth 



FIG. 6: Plot of the apparent VDF on a logarithmic scale for 
To = 30 eV, k = W and for different values of a = 5, 10, 20 
eV. Vth stands for the thermal velocity for T = 30 eV. These 
plot show the asymptotic power law behavior. The value of 
the exponent is -3.2 here. The dotted line corresponds to the 
Gaussian Doppler profile obtained for 30 eV. The deviations 
from this gaussian profile becomes more and more important 
as a grows. 
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FIG. 7: Plot of the apparent VDF on a logarithmic scale 
for To = 30 eV, ct = 10 eV and for different values of k, 
k — 1, 1/2, 10. Vth again stands for the thermal velocity for 
T = 30 eV, and the dotted line corresponds to the Gaussian 
Doppler profile obtained for 30 eV. 

field. The Fourier transform fa{k) of the apparent VDF 
is given by 

fa{k)^£ £„,_i(r) exp(^-^fc2^ dT, (40) 
and can be calculated explicitly using the following result 



[32^ which gives the Laplace transform of a Levy distri- 
bution 

p + co 

/ Cc,-i(T)exp{-sT)dT ^exp-cs", (41) 

with s = /2m + tO, m standing for the emitters mass. 
The apparent VDF is thus found to be a symmetrical 
Levy distribution of indexes a' = 2a and /?' = 

/^(-) = ^^-^0 . (42) 

Asymptotically, 

U^)-J^^ (43) 

in accordance with Eq. (|38|l . Here, the value of a is 
such that 1 < 2q; -I- 1 < 3 and therefore spans a different 
range than in the Sinai model. 

In this idealized model where only temperature fluctu- 
ates, the analysis of the apparent VDF tails, i.e. of the 
line wings, allows to retrieve information on the statisti- 
cal properties of temperature fluctuations. Indeed, power 
law decaying tails would for instance be a signature of a 
similar behavior for the temperature PDF. In addition, in 
this case, an analysis of the experimental value of the ex- 
ponent would allow to distinguish between different tur- 
bulence models, corresponding for example either to a 
Sinai or a Levy PDF. 

D. Discussion 

The study of the case where only one variable fluc- 
tuates leads to several enlightening conclusions. First 
of all, the Doppler profile is only affected by ion tem- 
perature and fluid velocity fluctuations along the line of 
sight, in contrast to the line brigthness which essentially 
reflects the variations of the density. In addition, the 
bulk of the line appears to be weakly sensitive to the 
presence of low frequency turbulence, unless the velocity 
fluctuations variance becomes comparable to the thermal 
velocity. Therefore, turbulence can indeed be neglected if 
we restrict ourselves to the study of the core of the line, 
as is usually done ITol . Conversely, the Doppler 
line wings behavior is significantly altered by turbulent 
fluctuations having non-Gaussian PDF. More precisely, 
long tails for the PDF translates into long tails for the 
apparent VDF, i.e. slowly decreasing line wings. In this 
sense, modifications on line wings are associated to in- 
termittency. A comparison with experimental spectra 
would require further work both from the theoretical and 
experimental sides, and will not be attempted here. In 
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particular, a refined model should simultaneously retain 
velocity and temperature fluctuations. Indeed, velocity 
and temperature effects cannot be distinguished a priori. 
In fact, examples where fluctuations of both fields lead to 
a power law behavior for line wings have been presented 
above. In addition, the couplings between density, veloc- 
ity and temperature fluctuations, which are responsible 
for anomalous transport, should also be taken into ac- 
count. To include these effects in our model, one could 
either rely on the determination of a joint PDF, or resort 
to a numerical integration of the fluid equations, which 
would allow a straightforward calculation of the appar- 
ent VDF from Eq. H14() . From the experimental point 
of view, line wings may seem difficult to measure, but it 
should be kept in mind that the acquisition time can in 
principle be chosen as large as needed. The only limita- 
tion here is the actual duration of the discharge station- 
ary phase during which the measurements are performed. 



be Gaussian. The main physical difference between the 
FPE and the FFFE given by Eq. iglj is the spatial non 
locality of the latter, obvious from the definition of the 
fractional derivative. This non locality is a consequence 
of the existence of flights connecting distant regions in the 
velocity space (the so-called Levy flights). This property 
can be traced back to the underlying description of the 
turbulent plasma. Indeed, at the microscopical scale the 
trajectory of the radiators can be modelled by a Langevin 
equation with gaussian white noise js^ . This model 
describes the collisional relaxation of the local velocity 
distribution toward the local Maxwellian Eq. {Tj). Using 
the fluctuation-dissipation theorem and the expression 
of the diffusion coefficient stemming from a random walk 
model ISa leads to 



(At;2 



kBT 



(46) 



VII. APPARENT NON-BOLTZMANN 
STATISTICS 

In the above section, we have shown that the apparent 
VDF may significantly differ from the Maxwellian calcu- 
lated using the averaged fields. For the sake of simplicity, 
let us only consider temperature fluctuations here. The 
fact that the apparent VDF can be a Levy distribution 
highlights a connection between spectroscopy, turbulence 
and anomalous statistics involving power-law tails, such 
as the Levy statistics. Indeed, it should be emphasized 
that in the case where no other observable than the spec- 
tral line shape is available (e.g. in Astrophysics), it is 
by no mean possible to determine whether the observed 
plasma is actually turbulent or homogenous. Therefore, 
if the temperature PDF is a Levy distribution Ca,~i{T), 
the Doppler spectra might be interpreted as resulting 
from an homogeneous and stationary plasma governed 
by Levy statistics. In other words, everything happens 
as if the plasma under study were in a non-equilibrium 
stationary state characterized by the Levy distribution 
of Eq. 1)42(1 . This stationary state can be seen as result- 
ing from a relaxation process governed by the following 
Fractional Fokker-Planck Equation (FFPE) [TllMlall 



dfa{v,t) - d 



D 



d\v\ 



2a ■ 



(44) 



where v and D are such that D/v — 2ac/{2m)". Here, 
the fractional derivative is defined in the sense of Riesz 



\krfeff 



(45) 



The usual Fokker-Planck equation (FPE) is recovered for 
a = 1. In our case a < 1, and the apparent VDF cannot 



where T is the local temperature and Tj the typical time 
between two jumps in the velocity space. Temperature 
thus determines the characteristic size of jumps in the ve- 
locity space. Therefore, high probabilities for large tem- 
perature fluctuations in the actual turbulent plasma im- 
ply high probabilities for flights in the apparent velocity 
space. This provides a simple physical picture explain- 
ing why the temperature PDF and the apparent VDF 
asymptotical behavior are linked, and leads to a deeper 
understanding of Eq. lp!H|l . Our results are reminiscent 
of those presented in references [s^ [s^, HI] , where a sim- 
ilar interpretation of Tsallis non extensive statistical me- 
chanics occurrence was proposed. The latter case arises 
if the ternperature PDF is such that 1/T is gamma dis- 
tributed |33- Let us emphasize that in our model, the 
temperature PDF shape is not arbitrary. In fact, it has 
to be determined from the fluid equation satisfied by the 
temperature field in the plasma under consideration, in 
which relevant expressions for both the source term and 
the fiux have to be specified (see Eq. O . For each of these 
expressions, the non-linear character of the latter equa- 
tion should give rise to a different non-gaussian statistical 
behavior, i.e. lead to a specific PDF, and therefore to a 
particular apparent statistics. A natural extension of this 
work would be to determine what properties fiuid equa- 
tions should have so as to lead to a Levy distribution for 
temperature. 



VIII. CONCLUSION AND PERSPECTIVES 

In this paper, we have presented a model retaining low 
frequency turbulence in Doppler line shape calculations. 
This approach is in particular relevant to the modelling 
of lines routinely measured in edge plasmas of fusion de- 
vices. We have shown that in presence of low frequency 
turbulence, a straightforward analysis of Doppler profiles 
yields an apparent velocity distribution function. This 
apparent VDF is a spatial and time average of the local 
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VDF. To investigate its shape, we have used a statisti- 
cal description of the plasma turbulent fluctuations, rel- 
evant whenever the acquisition time of the spectrometer 
is large with respect to the typical turbulent time scale. 
The resulting expression for the apparent VDF involves 
the joint Probability Density Function of the fluctuating 
fields. Next, considering the case where only one variable 
fluctuates, we have obtained several new results. While 
density fluctuations do not affect Doppler line shapes, 
velocity or ion temperature fluctuations can strongly in- 
fluence line wings. This is especially the case when their 
PDF have long tails such as power laws. It might there- 
fore be possible to diagnose such a behavior by the mean 
of line shapes, once Stark effect has been carefully ruled 
out. A reliable comparison with experiments would im- 
ply dedicated measurements which are not yet available, 
but also further modelling. In particular, the use of a 
turbulence code would be very helpful for diagnosis pur- 
poses, and this possibility will be investigated in a future 
work. Prom a more fundamental point of view, our work 



sheds light on some possible connections between turbu- 
lence, spectroscopy and non Boltzmann statistics, such 
as those involving Levy or Tsallis distributions. Our ap- 
proach furthermore relates the occurrence of one of these 
particular statistics to the properties of the fluid equa- 
tions describing turbulence. Our model thus provides a 
frame to investigate both experimentally and theoreti- 
cally some of the fundamentals aspects of the statistical 
properties of the physical observables in out of equilib- 
rium plasmas. 
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